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Abstract 



We discuss the collective modes in a harmonically trapped, highly- 
elongated Bose condensed gas. The transverse breathing mode exhibits a 
number of interesting features, such as the insensitivity of the condensate 
mode frequency to the interaction strength, and the closeness of the frequency 
to that of the non-condensed thermal cloud in the collisionless limit. Using 
finite temperature simulations, we show that these features are responsible 
for the very small damping rate observed experimentally. Our results for the 
temperature dependence of the damping rate and frequency shift are in ex- 
cellent agreement with experiment. We also demonstrate that the unusually 
small damping rate does not arise for the m = 2 mode or for more isotropic 
trap potentials, suggesting further possible experimental tests of our theory. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 
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I. INTRODUCTION 



The study of the finite temperature properties of trapped Bose-Einstein condensed gases 
touches on fundamental issues such as condensate growth and relaxation processes in an 
isolated quantum degenerate system. As such, they have been the subject of ongoing research 
since soon after the Nobel prize-winning discoveries of 1995. For example, early experiments 
at JILA Jl| and MIT || measured the frequency and damping rates of collective excitations 
as a function of temperature. More recent experiments have studied the finite temperature 
behavior of the scissors mode 0, and the radial breathing mode in a cigar-shaped condensate 
H. These experiments have in turn stimulated a considerable amount of theoretical work. 

Finite temperature theories have taken various forms. Early work focussed on solving 
Hartree-Fock-Bogoliubov (HFB) equations to find the excitation spectrum of the system 
0-^]. These excited states, when thermally occupied, determine the equilibrium density 
distribution of the non-condensed thermal cloud. Unfortunately, this method is unsatisfac- 
tory when considering the low energy collective modes since it treats the thermal cloud as 
static. In reality the condensate oscillations will induce fluctuations in the thermal compo- 
nent, which then act back on the condensate leading to damping and associated frequency 
shifts. One possible way of addressing this problem is to treat the thermal cloud dynamics 
perturbatively. This method has been pursued by a number of authors [PHT2"||, and although 
good agreement has been found with some of the experimental results , it is still unclear 
whether the method is useful in situations where the full dynamics of the thermal cloud is 
important. 

An alternative approach is to formulate a quantum kinetic theory, which in principle 
treats the dynamics of the system fully and consistently. This includes second-order colli- 
sional processes, which lead to equilibration of the system. Kinetic theories for dilute Bose 



gases have been formulated by Stoof ||13|| , Gardiner and co-workers ||14| , Walser et ai, (15 
and Zaremba, Nikuni and Griffin (ZNG) ||16|| . Finally, there is the possibility of treating the 
evolution of the system using Fokker-Planck equations for quasiprobability distributions. A 
specific example of this is the truncated Wigner representation, or classical field method 
[TT| , [TE| ], which is valid when the occupation of the relevant states is much greater than unity. 

In this paper we shall present results of solving the ZNG equations using a method 
detailed in a previous publication [jEJ]. Earlier applications of the method yielded results 
in good agreement with the scissors mode experiment at Oxford PJ2^| as well as with 
the quadrupole modes in the JILA experiment Here we focus instead on the ENS 

experiment [|J, adding more details to the discussion of our simulations in Ref. [22]|. The 
ENS experiment studied the radial breathing mode of an elongated condensate, and found 
unusually small damping rates and frequency shifts as compared to other experiments. We 
explain this behaviour in terms of an accidental degeneracy between the condensate and 
thermal cloud oscillation frequencies. In Section II we review the ZNG theory and our 
numerical scheme, while in Section III we discuss the collective modes in an elongated 
harmonic trap. In Section IV we present results from our finite temperature simulations 
for the tranverse breathing and quadrupole modes in this geometry. For the latter, the 
condensate and thermal cloud frequencies are no longer degenerate and we find a much 
larger damping rate and frequency shift. We also show that larger damping rates arise when 
the degeneracy is lifted in more isotropic trap potentials. Finally, Section V presents our 



2 



conclusions. 



II. ZNG THEORY 

We first give an overview of the ZNG theory which forms the basis of the numerical 
simulations described in this paper. Further details of the formalism, including the ap- 



proximations made, can be found in Ref. [16]. In the theory two coupled equations are 



used to describe the condensate and thermal cloud dynamics. In particular, the condensate 
wavefunction $(r, t) evolves according to the generalized Gross-Pitaevskii (GP) equation 

d ( h 2 V 2 \ 

iK— $(r,t)= ( + U ext {r,t)+g[n c {r,t) + 2h{r,t)] -iR(r,t)J $(r,t). (1) 

The condensate interacts locally with other condensate atoms and with the non-condensate, 
with interactions parametrized by g = 4irh 2 a/m, where a is the s-wave scattering length 
and m is the mass of the atoms. These mean-field interactions involve the condensate 
density n c (r,i) = |$(r, t)| 2 and the non-condensate density n(r,t). The term U cxt (r,t) = 
m[ujj_(x 2 + y 2 ) + uj 2 z 2 \l1 represents the trap potential. The frequencies in the radial (lj±) 
and axial (u z ) directions are generally different, and the ratio A = lu z /u± conveniently 
parameterizes the trap anisotropy. 

The thermal cloud is described in terms of a semiclassical Wigner distribution /(p, r, t) 
which defines the non-condensate density according to 

S(r - i) = /(2^ /(p ' r - t) - (2) 
This distribution function satisfies the kinetic equation 

|-/(p, r, t) + P ■ V/(p, r, t) - W(r, t) • V p /(p, r, t) = C 12 [/] + C 22 [/] , (3) 
at m 

where the term U(r,t) = U cxt (r) + 2g[n c (r,t) + n(r, t)] reflects the combined effect of the 
trap and mean-field potentials on the thermal cloud evolution. Interactions also enter in 
the terms on the right-hand side, which are of order g 2 and correspond to binary collision 
events between thermal atoms (C22) and between thermal and condensate atoms (C12). The 
second process is particulary important in the context of condensate growth, since it leads 
to a transfer of atoms between the condensate and thermal cloud. It also accounts for the 
non-Hermitian term in the GP equation (Q) which is given by 

This term has the effect of changing the normalization of the condensate wavefunction. 

We note that the kinetic equation (^) for the thermal excitations is the same as that 
for gas of classical particles moving in the effective potential U. This is the basis of our 
numerical solution of this equation, which tracks the Hamiltonian dynamics of a swarm 
of test particles using a sympectic integrator. Collisions are treated using a Monte Carlo 
evaluation of the collision integrals. The GP equation can be propagated in time using a 
variety of methods, but we favour a Fast Fourier Transform split-step operator technique. 
More details of the numerical methods can be found in Ref. [19[ . 
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III. COLLECTIVE MODES 



In this section we discuss the collective modes of the system, with particular emphasis 
on the so-called transverse breathing mode in a cigar-shaped trap. We begin by reviewing 
results for a pure condensate at temperature T = 0. In this limit the condensate dynamics 
is described by (]l|) with h = and R = 0, i.e., 



ih— $(r,t) = l-—— + U cxt (r)+gn c (r,t)\ $( 



,t). (5) 



The equilibrium solution can be found from (|5|) by making the substitution $(r, t) = 
^o(r)e~ tflt , where /i is the chemical potential. When the number of atoms is sufficiently 
large, the kinetic energy term is dominated by the other terms and can be neglected. This 
defines the Thomas-Fermi (TF) approximation and gives an equilibrium density profile of 

n c0 (r) = - \p - C/ ext (r)] . (6) 

The frequencies of condensate excitations can be found by linearizing (|), and analytic 



solutions for the lowest modes in the TF limit were first presented by Stringari ||23|| . In a 
spherical trap (A = 1), the modes can be labelled by the angular momentum quantum num- 
bers (l,m). The lowest of these modes correspond to monopole (0,0) and quadrupole (2,0) 
excitations with frequencies ^/5to± and \[2{ J j± 1 respectively. In an anisotropic, axisymmetric 
trap these two modes are coupled to give the dispersion relation 
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w| = u)\ \ 2 + -A 2 ± -V9A 4 - 16A 2 + 16J . (7) 

In Fig. 1 we plot these two solutions as a function of A. The two curves tend to the monopole 
and quadrupole frequencies of the isotropic trap in the A — > 1 limit. For the cigar-shaped 
condensate which will be of interest here (A < 1), the low frequency mode u>- corresponds to 
a breathing-like oscillation predominately along the axial (z) direction, while the high-lying 
mode u + is associated with the radial breathing mode. In the limit A — > the frequency of 
the latter approaches 2u±. For example, in the experiment of Ref. 0, A ~ 6.46 x 10 -2 and 
the transverse breathing mode has a frequency of uj + ~ 2.00052u;_|_. 

These results are valid in the strongly interacting TF regime, but the non-interacting 
limit g = is also instructive. In this case the transverse quadrupole and breathing modes 
mentioned above both have a frequency of 2u>j_. This clearly has interesting consequences 
for the breathing mode in the A = limit, since it implies that the frequency is the same 



in the non-interacting and strongly-interacting limits. Indeed, calculations confirm [24| that 
this is also true in the intermediate regime, demonstrating the unusual property that the 
frequency of this mode is independent of interactions. This behavior is reminiscent of that 



for a gas in a 2D harmonic trap interacting by means of a contact potential |p5jp6| , where 



there exists an undamped breathing mode with frequency 2u±, independent of interaction 
strength, temperature, or statistics. 

The thermal cloud dynamics above T c is given by (0) with n c = and Cyi = 0, i.e., 

-f(p, r, t) + ^ • V/(p, r, t) - W ext (r, t) ■ V p /(p, r, t) = C 22 [f] , (8) 
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where for simplicity we have also neglected the thermal cloud mean field term, 2gn, since it 
is has a very small effect on the dynamics in the normal state. The lowest collective modes 
for the thermal gas can be conveniently found from (El) by taking moments. This procedure 
was used to study a classical gas in p7| , though the results for a non-condensed Bose gas 
are essentially the same. The frequency and damping of the modes depend upon the mean 
collision rate r _1 relative to the frequency of the mode u. In the purely collisionless regime 
lot — > oo, and the frequencies of the transverse quadrupole and breathing modes are both 
to = 2uj±, while the damping is zero. In the hydrodynamic limit (ut — > 0) the modes 
are also undamped, but here the frequencies of the respective modes are luq = ^/2cuj_ and 
lu m = a/10/3u;_i_ (where the latter is specifically for the cigar-shaped limit A <C 1) p7|j28|| . In 
general the modes will be damped through collisions, and will have a frequency intermediate 
between these two limits. In the experiment of [|J the system resided in the near-collisionless 
regime, so that both thermal cloud modes were only weakly damped with frequency uj ~ 2u±. 

In order to understand the transverse breathing mode more fully it is instructive to derive 
an equation of motion for the operator R 2 = Yui( x i + vf)- The expectation value of this 
operator for a particular dynamical state gives the mean-squared radial size of the cloud 
and it is therefore a suitable variable to describe the transverse breathing mode of interest. 
For the Hamiltonian 



we have the Heisenberg equation of motion 



^rr^ Q 



, -r[H,R*\ = 2-, (10) 
dt n m 

where the operator Q is given by 

Q = y^jXiPix + ViViy + h.C.) . (11) 

i 

After differentiating ([TOD again and evaluating the commutator [H, Q] we obtain the result 

^ + (2u ± ) 2 R 2 = ^(H — T z — U?% (12) 

where T z = ^2i{p 2 z /2m) is the kinetic energy associated with motion in the ^-direction and 
U z xt = J2i(, mu}2z i It is worth emphasizing that this equation for R 2 is an operator 
identity which is equally valid classically as it is quantum mechanically. 



We now consider the expectation value of (0) with respect to some nonequilibrium 
density matrix. Denoting expectation values by angular brackets, and noting that the energy, 
E, of the system is a constant of the motion, we have 

^ + {2u i _) 2 x = ~[E - (T z ) - (Ur)] , (13) 
dt z m 

where we have defined \ = (R 2 )- For a strictly two-dimensional (2D) system, (T z ) = 
(U* xt ) = 0, and (0) has the solution 
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X = Acos(2u ± t) + x- (14) 

In other words, x oscillates about x — E/(muj'j_) with frequency 2u±. This behaviour was 
first noted by Pitaevskii and Rosch in their analysis of the 2D system. For three-dimensional 
systems, (T z ) and (U z xt ) are themselves dynamical variables. Their appearance in ([13]) will 
in general lead to a shift of the frequency of the transverse breathing mode from 2lj±, except 
in certain limiting cases. 

The equation of motion for mean values in (|T3|) is also valid when the dynamics is 
governed by the time-dependent GP equation. However in this case, the average of an 
operator O is defined as (O) = J dr $*(r, t)OQ(r, t). In the limit of a cylindrical trap, 
A — > and (U z xt ) = 0. If we then consider an axially symmetric mode which is spatially 
invariant along the length of the trap, we have (T z ) = and ( |i~3"D implies that a mode at 
the frequency 2u± exists. This particular mode was previously identified as the transverse 
breathing mode in a cylindrical trap. It should be emphasized that this conclusion is based 
on the approximate GP description. In the many-body treatment, (T z ) ^ and the existence 
of a 2u± mode then depends on (T z ) being constant. It is not known whether a dynamical 
state with this property exists in general. 

In applying ( |T3| ) to a semiclassical gas we interpret the expectation value of a physical 
observable O as the phase space average (O) = J(dpdr//?, 3 )(9(p, r)/(p, r, t). Two limits 
are of interest: the collisionless and the hydrodynamic. In the former case one can use (^) 
to show that d((T 2 ) + (U z xt ))/dt = for a separable confining potential. Eq. (]i~3|) then 
implies that a collisionless gas has a transverse breathing mode at the frequency 2u±, as 
noted earlier. 

Once collisions are included, ((T z ) + (U z xt )) is no longer a constant of the motion and 
the frequency of the transverse breathing mode will change. We can illustrate this for a 
cylindrical trap (u z = 0) in the hydrodynamic limit. In this regime, collisions are sufficiently 
frequent to ensure local thermodynamic equilibrium, so that the distribution function is of 
the Bose form 

/(P) r > t) = e ^( rit ){[p_ mv ( r|t )]2/2 m - (1 (r | t)}-l ' (15) 

where /3(r,t), /x(r, t) and v(r, t) are all local hydrodynamic variables. For the transverse 
breathing mode the local velocity field has the form v(r, t) = ar±, where r_|_ is a vector in the 
radial direction, tranverse to the axis of the trap. Using (|i~5|) one can show that (T z ) = (T)/3 
and E = (T) +ma;^x/2, where we have neglected the nonlinear contribution of the velocity 
field to the kinetic energy (T) = (p 2 /2m) of the mode. By substituting these expressions 
into ( |T3"D one finds that the transverse breathing mode has frequency u>m — a/10/3o;_i_, m 



agreement with other calculations in this limit p7| , |28[| . This provides a concrete example of 
a situation in which a non-constant (T z ) leads to a shift in the mode frequency from 2uj±. 
It is therefore clear that a mode at this frequency is not a universal property of a trapped 
gas in the cylindrical geometry. 

IV. NUMERICAL RESULTS 

The behavior of the condensate at < T < T c is influenced by the presence of thermal 
excitations. Normally this leads to a decay of condensate oscillations, which in the near- 
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collisionless regime arises through Landau damping. There is also an associated frequency 
shift. However, in the experiment of the transverse breathing mode was found to be 
very weakly damped, with a frequency almost independent of temperature. This is in stark 
contrast to the much larger damping rate found in their measurements of the quadrupole 
mode, as well as in other experiments Jlj-Q. In this section we present results of both 
the breathing and quadrupole modes for the system studied in Ref. ||. We show that the 
observed behavior of the breathing mode is a consequence of the degeneracy of the condensate 
and thermal cloud mode frequencies coupled with the fact that both components are initially 
excited. 

The experimental scheme used to excite the system was to suddenly change the radial 
trap frequency and then to reset it to its original value after some short time r. This stepped 
excitation can be represented by 

cu' x (t)=uj{l + a[e(t)-e(t-T)}}, (16) 

where we take u± = 2tt x 182.6 Hz, a = 0.26, and uj±t = 0.172. Although this is a relatively 
simple method for exciting condensate oscillations, in finite temperature studies it has the 
disadvantage of also exciting the thermal cloud. This aspect is relatively unimportant when 
the oscillation frequencies of the condensate and thermal cloud are sufficiently different, but 
when they are similar it can have a profound effect on the behavior of the system. An 
important implication is that the usual perturbative calculation of Landau damping which 
assumes the thermal cloud to be at rest [§-[l2| are no longer applicable. 

To illustrate this point, we show in Fig. |2|(a) the result of a simulation at T = 125 nK (as 
compared to the experimentally measured critical temperature of T c ~ 290 nK). We plot the 
radial moments for the condensate, Xc, and thermal cloud, Xn, as a function of time, where 
each has been divided by its initial value to give the relative oscillation amplitude. One 
sees that both components respond to the stepped excitation by oscillating in phase with 
approximately equal amplitudes. In addition, the oscillation frequency is very close to 2u± 
and the damping rate is very small. As discussed earlier, the frequency is that expected for 
the condensate for A ~ and for the thermal cloud in the collisionless regime. The common 
in-phase oscillation of the two components is further reinforced by mean-field interactions. 
The lack of relative motion between the two components accounts for the fact that Landau 
damping is not effective. Evidence supporting this interpretation is provided in Fig. @(b), 
where we plot the result of exciting only the condensate by imposing a radial velocity field, 
v c oc xi + yy In this case the condensate oscillates initially in the presence of a stationary 
thermal cloud and as a result, Landau damping is fully active, giving rise to the much larger 
damping rate seen here. 

To quantify the frequency and damping rate, we fit the condensate oscillation to an 
exponentially decaying sinusoid Acos{uit + <p)e~ Tt + C. Fig. |3| compares our results at 
different temperatures with the experimental data. We see very good agreement between 



simulations and experiment for both the frequency and damping rate. As discussed in jzz 
the damping that occurs is essentially due to Cyi and C22 collisional processes. This is unlike 
the situation for other modes in the near-collisionless regime |I^.2l]fl where in fact Landau 



damping is the dominant damping mechanism. As indicated in Fig. 0(b), Landau damping 
of the breathing mode only appears when the two components are moving relative to each 
other. 
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Recently, results from a calculation by Guilleumas and Pitaevskii appeared [[29] that 
point to a different conclusion. They performed a perturbation theory calculation similar to 
|TT| , but for a cylindrically trapped condensate. In this geometry the excitation spectrum, as 



obtained from the Bogoliubov equations, is characterized by a series of curves continuous in 
k, the axial wavevector, with each curve corresponding to distinct radial (n) and azimuthal 
angular momentum (m) quantum numbers ||30|| . Calculation of the Landau damping rate 
then consists of evaluating the following expression 

r = ^2 IA?| 2 (/< ~ fj^if^J ^osc) , (17) 

ij 

where the delta function reflects energy conservation of the decay process, in which a quan- 
tum of oscillation ftw osc is annihilated, and the z-th excitation transformed into the j-th. 
The probability of this transition is determined by the matrix element (as defined in 
Ref. [^9|), and the distribution function /. The thermal cloud is assumed to be in thermo- 



dynamic equilibrium, so that the states i, j are thermally occupied according to the Bose 
distribution fi = [exp(ftwi/kBT) — l] -1 . Using this method they calculated damping rates 
for the lowest m = and m = 2 modes in the k = limit, which corresponds to the trans- 
verse modes in the elongated 3D geometry. We shall discuss the m = 2 results later, but it 
is interesting that they found a damping rate for the m = which is an order of magnitude 
smaller than the experimentally measured value. Since this calculation assumes a stationary 
thermal cloud, it is in fact two orders of magnitude smaller than our result for the analogous 
situation where only the condensate is excited initially [Fig. @(b)]. 

The discrepancy between the two results could arise for a variety of reasons. It may 
be that the semiclassical approximation, as employed in our calculations, is not valid for 
this particular case. The small damping rate in Ref. |[{| is partly the result of the small 



number of transitions that satisfy energy conservation. This constraint is less severe in the 
semiclassical approximation in which the excitations are continuous in both the axial and 
radial directions. The damping rate seen in [Fig. 0(b)] might therefore be an overestimate. 
On the other hand, there is also the possibility that the small damping rate is an artifact of 
the calculational method employed for the infinite cylinder. In the spherical trap discussed 
by the same authors |Tl[ , energy can never be conserved precisely in a transition due to the 



discreteness of the energy levels. So in practice the Landau damping rate (\L7\) is calculated 
by replacing the delta functions with Lorentzians of finite width. In a previous publication 
|T9| we found good agreement between this approach and our simulations. This procedure 
would have to be employed for any 3D system with a finite value of u z . Since the excitation 
spectrum is a continuous function of k for an infinite cylinder, this procedure for dealing 
with discrete excitations is not required. The neglect of off-energy shell transitions in this 
case may be leading to an underestimate of the damping rate as compared to the finite trap 
situation. In order to address these questions it would be useful to extend the perturbative 



calculations of [11] to elongated condensates. One could then study the variation of the 
damping with anisotropy to see whether the results begin to diverge from the semiclassical 
simulations at some point. We are currently exploring this possibility. 

As we discussed above, the small damping rates and frequency shifts shown in Fig. |3] are a 
consequence of the fact that there is no relative motion of the condensate and thermal cloud 
in the transverse breathing mode. This behaviour does not occur if the two components 



S 



tend to oscillate at different frequencies as is the case for the m = 2 transverse quadrupole 
mode. To demonstrate this we have performed simulations for the m = 2 mode by imposing 
identical velocity fields of the form V{ c ,n} oc xi — ?/j on the two components. An example of 
the subsequent oscillations of the condensate and thermal cloud is shown in Fig. |j. One sees 
that the frequencies are indeed quite different in this case - the condensate oscillates with a 
frequency approximately equal to v^2^±, while the thermal cloud oscillates at a frequency 
close to 2uj^. As for the m = mode, we fit the condensate data to a damped sinusoid, and 
plot the frequency and damping rate as a function of temperature in Fig. |5|. In contrast to 
the m = data in Fig. the m = 2 mode shows a much larger shift in frequency, as well as 
an order of magnitude increase in damping rate. The relative motion between the condensate 
and thermal cloud in this case accounts for these increases. Indeed, from the perspective 
of the condensate, the thermal cloud motion averages over the course of the simulation to 
that of a stationary cloud. This is confirmed by simulations where the thermal cloud is not 
excited initially. As seen in Fig. [|, these simulations exhibit very similar results to when 
both components are perturbed, especially when T < 150 nK. The differences between the 
two schemes at higher temperatures mainly reflect statistical variations of the results from 
one simulation to the next. 



Interestingly, when Guilleumas and Pitaevskii [29] investigated the m = 2 mode for the 
cylindrical limit, they found damping rates in line with our results in Fig. [5], especially at 
higher temperatures. There are some discrepancies at low temperatures, but it must be 
remembered that precise agreement would not be expected between the two calculations 
due to differences in geometry, parameters and calculational method. Nonetheless, it is 
intriguing that such good agreement should be found here and not for the m = mode. 
Again, further study is needed. Experiments for the m = 2 mode would also be of interest, 
since comparison to the results in Fig. |5] would provide further information regarding the 
reliability of our theory. 

Finally, we discuss another method for breaking the degeneracy of the condensate and 
thermal cloud modes. As seen in Fig. [I], increasing the parameter A leads to a rise in 
the breathing mode frequency u + , while at the same time the thermal cloud frequency 
remains near 2u±. Thus, as the trap becomes more isotropic the condensate frequency 
moves away from degeneracy, and one would expect the damping of the breathing mode to 
correspondingly increase. To check this we have performed simulations for larger values of 
A. An example of the response to a stepped excitation fll6|) is shown in Fig. |] for A = 0.75; 
the geometric mean of the trap frequency is chosen to be the same as in the ENS experiment 
(u = \ 1 / 3 uj± = 2n x 73.3 Hz). Fig. ||(b) shows the time dependence of the thermal cloud 
moment Xn- It exhibits a damped oscillation at essentially one frequency; a single- mode fit 
to the data yielded a frequency of u/u± ~ 2.007 and a damping rate of T/u± ~ 0.019. The 
latter is approximately four times larger than found for the ENS geometry. The behavior 
of Xc in Fig- 0(a) is more complex due to the fact that both of the u± modes are being 
excited. In fact, a combination of two damped sinusoids provided a very good fit to the time 
dependence shown and yielded the following fit parameters: u_/u± ~ 1.12, T_/uj± ~ 0.037; 
u+/uj± ~ 2.03, T + /u± ~ 0.024. The frequencies are quite close to the values given by (^) 
for A = 0.75. More importantly, we see that the relative damping rate of the u + mode is 
approximately an order of magnitude larger than that of the transverse breathing mode from 
which it evolves. Thus even a relatively small difference in the frequencies of the condensate 
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and thermal cloud modes is sufficient to significantly enhance the Landau damping rate. 
This fact emphasizes the very special nature of the transverse breathing mode in the A — > 
limit. 

V. SUMMARY 

In this paper we have investigated the transverse oscillations of a Bose gas in the context 
of semiclassical finite temperature theory. Analytical results were derived in the cylindrical 
limit for the dynamics of the condensate at T = and thermal cloud for T > T c . Simulations 
are used to study the dynamics in an elongated trap at < T < T c , for the parameters 
of the experiment of ||. Our results for the frequency and damping of the m = mode 
demonstrate very good agreement with this experiment. However, comparison to the results 
of a related theoretical study in |29| raises several interesting questions. We also study the 
m = 2 mode, and find a much larger damping rate and frequency shift associated with the 
presence in this case of relative motion between the condensate and thermal cloud. Finally, 
we study the effect on the m = mode of decreasing the trap anisotropy, and find that the 
damping increases dramatically as soon as the accidental degeneracy between the condensate 
and thermal cloud frequencies is lifted. 
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FIG. 1. Frequencies of the lowest m = condensate modes in the Thomas- Fermi limit, 
given by (0), as a function of the anisotropy parameter A = lo z /uj_\_. The higher mode u;+ (which 
corresponds to the transverse breathing mode in the A — > limit) is labelled by the solid line, while 
uj- is denoted with the dashed line. 
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FIG. 2. Time-dependent radial moments for the condensate, Xc (narrow line) and thermal 
cloud, Xn (bold line), divided by the corresponding values at t = 0. The figures are for a tempera- 
ture of T = 125 nK, and show the result of (a) exciting the system using the "tophat" perturbation 
scheme employed experimentally ([R]), and (b) exciting the condensate only by initially imposing 
a velocity field of the form v c oc x'\ + yj . 
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FIG. 3. Frequency (a) and damping rate (b) of the condensate transverse breathing mode. 
Our results (open circles) are compared to the experimental data of [|j (solid circles), where the 
simulation parameters and excitation scheme are chosen to match the conditions of the experiment. 
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FIG. 4. Time-dependent quadrupole moments (x 2 — y 2 ) for the condensate (bold line) and 
thermal cloud (narrow line), where the moments have been scaled so that they overlay one another. 
The plots are for a temperature of T = 125 nK, and show the result of initially imposing a velocity 
field of the form V{ c n } oc xi — yj on both components. 
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FIG. 5. Frequency (a) and damping rate (b) of the condensate transverse quadrupole mode, 
where we show the results of initially exciting both components (open circles) and the condensate 
only (bullets) by imposing a velocity field. Simulation parameters are the same as in Fig. ||. 
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FIG. 6. Time-dependent radial moments for (a) the condensate and (b) the thermal cloud, 
divided by the corresponding values at t = 0, for a less anisotropic trap. The figures are for 
a temperature of T = 125 nK, and show the result of perturbing the system using the stepped 
excitation in (|16|). The trap parameters are A = 0.75 and ui = 2ir x 73.3 Hz. 
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